Helical liquids and Majorana bound states in quantum wires 
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We show that the combination of spin-orbit coupling with a Zeeman field or strong interactions 
may lead to the formation of a helical liquid in single-channel quantum wires. In a helical liquid, 
electrons with opposite velocities have opposite spin precessions. We argue that zero-energy Ma- 
jorana bound states are formed in various situations when the wire is situated in proximity to a 
conventional s-wave superconductor. This occurs when the external magnetic field, the supercon- 
ducting gap, or the chemical potential vary along the wire. We discuss experimental consequences 
of the formation of the helical liquid and the Majorana bound states. 



States of matter that support Majorana fcrmions have 
received much attention in the context of quantum com- 
putation. A widely separated pair of Majorana bound 
states forms a nonlocal fcrmionic state which is immune 
to local sources of decoherence, thus providing a platform 
for fault-tolerant quantum memory. Moreover, since Ma- 
jorana states realize a representation of the non-Abelian 
braid group, topological quantum information processing 
can, in principle, be effected by braiding operations 
A realization of such states where they can be readily 
moved around and manipulated is therefore highly desir- 
able. 

There are several suggestions for physical systems that 
support Majorana states, for particular ways to manip- 
ulate them, and for measuring their properties. These 
include fractional quantum Hall states at filling factor 
v = 5/2 Q, p-wave superconductors 0|, surfaces of 
three-dimensional topological insulators in proximity to 
a superconductor [4(, and helical edge modes of two- 
dimensional topological insulators in proximity to a fer- 
romagnet and a superconductor [f|. Recently, it was 
suggested that a semiconducting thin film sandwiched 
between an s-wave superconductor and a magnetic insu- 
lator Q will host Majorana states associated with super- 
conducting vortices. All these proposals, however, would 
be experimentally extremely challenging. 

Realizing and manipulating Majorana fermions in 
wires may be decisively simpler. Here we show that quan- 
tum wires with strong spin-orbit coupling, such as InAs 
or InSb wires, and banded carbon nanotubes form a he- 
lical liquid, similar to the edges of a topological insula- 
tor. Consequently, these wires support Majorana states 
when they are in proximity to s-wave superconductors, 
and a magnetic field. Most importantly, we explain how 
they can be produced and manipulated by variations of 
a chemical potential, which could be simply produced 
by a set of micron-sized gates capacitatively coupled to 
the wire. Below we outline the key physical properties 
of Majorana states in quantum wires, their experimental 
signatures, and how Majorana-supporting wires could be 
extended into networks of Majorana fermions, enabling 



quantum information processing. 

Our analysis begins with writing the Hamiltonian for 
a spin-orbit coupled quantum wire. Without loss of gen- 
erality, let us choose the wire to lie along the y direction, 
the spin-orbit interaction, u, to be along the z direc- 
tion, and a magnetic field B to be along the x direction. 
In addition, the wire is in contact with a superconduc- 
tor, with the proximity strength being A (assumed to be 
real). The Hamiltonian is given by [J] 

H= [&(y)H*(y)dy; * f = (1) 
H = [p 2 /2m - fi(y)] t z + u(y)p <j z t z + B(y)a x + A(y)r x . 

Here, V , t,(4-)(f) ^ s ^ ne annihilation operator of electrons 
with up (down) spin at position y. The Pauli matrices 
a and r operate in spin and particle-hole space, respec- 
tively, fi is the chemical potential. 

In the absence of the Zeeman field and the supercon- 
ducting proximity, the eigenstates of the Hamiltonian (JTJ 
will have an energy-momentum dispersion consisting of 
two shifted parabolas crossing at momentum p = 0. The 
Zeeman field B removes the level crossing and opens a 
gap at p = 0. We note that such a gap may also occur 
due to strong electron-electron interactions @, H[, and 
therefore B should be generally contrued as either a mag- 
netic field perpendicular to the spin orbit coupling, or an 
interaction induced gap. The pairing A will play two 
crucial roles: Opening a gap at the outer wings of the 
dispersion, where the Zeeman field is unimportant, and 
modifying the gap forming near p = 0. The former role 
eliminates the possibility of high-momentum gapless ex- 
citations, thus leaving only the chiral states near p = 
as low energy excitation. These states resemble the edge 
of a topological insulator 0, The latter role allows 
us to tune the topological phase transitions essential for 
the production of Majorana fcrmions. Note that another 
way for gapping out the large momentum excitations is 
by coupling our system to an antiferromagnet with peri- 
odicity comparable to 2k p of the wire. Interactions may 
also open a pairing ifc^gap for chemical potentials away 
from the Zeeman gap 
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The emerging spectrum for constant /i, u, A, and B, 
is conveniently obtained by squaring the Hamiltonian 
twice. This straightforwardly yields the expression: 

Eg = B 2 + A 2 + £ + (up) 2 ± 2^B 2 A 2 + B 2 e p + {up) 2 e p 

(2) 

where £ p = p 2 /2m — p. Fig. [TJ displays the spectrum for 
several values of £?, A, and p. As these parameters vary 
(while B and A remain nonzero), a gap closing and re- 
opening indicates a topological phase transition. Gener- 
ically, we expect gaps appearing near p = and near the 
Fermi momenta corresponding to £ p ± up = 0. We will 
denote these gaps as Eq and E\, respectively. 




FIG. 1. (a) Single-particle spectrum for u = B = A = 0. 
(The two colors denote the different spin components). The 
energy (momentum) scale is set by mu 2 /2 (by mu), with u 
the spin-orbit coupling strength, (b) Excitation spectrum of 
adding or removing an electron for n — B — A = 0. (c) 
Excitation spectrum for B — 1/4, A = p = where a spin 
gap opens near p = 0. (d) B = 2.5, A = 1/2, p, — with a 
superconducting gap in the wings and a spin gap near the 
origin. This situation is analogous to a p-wave supercon- 
ductor. We refer to this phase as the "spin gap phase" (e) 
B = 1/4 = A = 1/4, fi = 0. The gap near p = closes, 
the gap at finite p persists. At this critical point a quan- 
tum phase transition occurs, (f) B = 1/4, A = 0.3, p, — 0. 
All gaps in the excitation spectrum are controlled by A. (g) 
B = 1/4, A = A,a = VB 2 - A 2 = v2T/20. The gap at 
p = closes due to the shift in chemical potentia 1. (h) 
A superconducting gap opens up in the entire spectrum due 
to the shift of the chemical potential above its critical value 
A = 1/10, B — 1/4, n = 0.3. 



As hinted above, it is the zero- momentum gap, Eq, 
which is crucial for our understanding of the emerging 
Majorana states. Examining E- at p = we notice that 

E = E(p = 0) = \B-^A2 + ui\. (3) 

For B 2 > A 2 + p 2 , Eq is a spin gap due to the Zeeman 
field (or strong interaction), while for B 2 < A 2 + p 2 it 
is a superconducting gap, thus when B 2 = A 2 + p 2 a 
quantum phase transition occurs. At the same time the 
gap Ei near p 2 = 2pm is always a superconducting gap, 
as we require A to always remain finite. 

The phase transition evident in Eq allows the forma- 
tion of Majorana states. Indeed, the dependence of Eq 
on B, A, and p enables us to construct zero-energy Ma- 
jorana states in various ways. As in edge states of 2D 
topological insulators a Majorana bound state will 
form when B changes in space and crosses A, e.g. at 
y = (cf. Fig [5b), or when A varies in space and crosses 
B (cf. Fig [2H). 

Here wc emphasize, however, a third possibility: vary- 
ing the chemical potential, p. Let us assume that B > A 
so that for p = we have a spin gap Eq. But when 
p > \JB 2 — A 2 , the gap Eq, Eq. ©, is clearly supercon- 
ducting. Thus, we can form a Majorana state by tuning 
p between these two values (cf. Fig [2p) . We note that 
changes in p do not significantly influence the gap E\ , so 
that the electronic states near ±&f do not play a role. 

The one-dimensional geometry allows for a simple 
demonstration of how to form Majorana states where 
their wave functions can be obtained essentially exactly. 
Let us consider these examples in a long ring with one 
conducting channel, in proximity to a superconductor 
and a Zeeman field, as illustrated in Fig. [2^. Since the 
relevant momenta are near p = 0, in the treatment below 
we use the Hamiltonian linearized in that region: 

H = up a z r z - p{y)r z + B(y)a x + A(y)r x (4) 

Spatially varying B. Assume A > is constant, 
p = 0, and that B > A for y > and B < A for y < 
(Fig. [2Jd; note that the periodic boundary conditions re- 
quire another point where B = A). Near the crossing 
point y = 0, we write B(y) = A + by. Due to particle- 
hole symmetry, it is useful to square the Hamiltonian 
Eq. ((4]) to diagonalize it. In addition to the square of 
each term and the mixed BA term, we obtain a term 
{upo~ z T z , Bo~ x } = ia y T z u[p, B] = a y T z ub which arises be- 
cause B depends on space and does not anticommutc 
with the spin-orbit coupling. Collecting all terms, we 
have 

n 2 b = (up) 2 + B(y) 2 + A 2 + uba v r z + 2AB(y)a x r x (5) 

Rotating Hi by Uj = 1/2 (r 2 - ir x - ia x j z + a x r x ), we 
find that Ub ■ Hi ■ is diagonal with components 
(up) 2 + (A ± B) 2 ± ub. The interesting modes are those 
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with a minus sign in the brackets, A — B. They corre- 
spond to a simple harmonic oscillator Hamiltonian with 
ground-state wave function <p(y) = (6/(ii7r) 1 / 4 )e -&2 '"/( 2tl ' 

and energies E 2 = 2ub(n + 1/2) ± ub, n = 0,1,2, 

For b > 0, the minus sign yields a zero-energy state with 
Bogoliubov operator 

7b = ib = (m - m) = g {h - % + #1 + V't) 

f/i = 1/V2 (v4 + ^t) 5 % = l/(V2i) (Vl - ^) ■ (6) 

The Majorana state at the second crossing point along 
the ring follows by b — > —b. Thus, this zero-energy state 
is Eq = with Majorana operator — i/\/2(r?i + Tj2)- 

Spatially varying A. For the case where A depends 
on y, we assume A(y) = B + dy, (i = 0, and a constant 
B (Fig[2t). The Hamiltonian here is similar to that in 
the y-dependent B case, if we exchange r and a in Eqs. 
((H) and Therefore, the Majorana states emerge in 
this case in exactly the same way as above, except with 
the diagonalizing matrices being Jj\ = Ul{r <h> a), and 
with b and A exchanged with d and B respectively in 
the resulting wave function. This yields (for positive d) 

Id = 7d = (vi - m) A/2. 




B-A 



,=fi. 



FIG. 2. (a) Wire in a ring geometry. Both halves have con- 
stant parameters and are joined by short junctions with a 
linearly varying parameter. Majorana states (marked by cir- 
cles) are formed at the junctions, (b) Majorana state in the 
sector p = when B varies. The gap in the finite-p sector 
remains finite in the entire wire, (c) Majorana state in the 
sector p — when A varies, (d) Majorana state in the sector 
p = when p varies, (e) "p-wave" Majorana state when A 
changes sign. The sector p = remains gapped in the entire 
wire. Each crossing with A = hosts two Majorana states. 

Spatially varying /i. If B > A in the entire wire, 
then at the interface between spin-gap regions with /i 2 < 
B 2 — A 2 and pairing gap regions with \fi\ 2 > B 2 — A 2 , a 
Majorana state will also form (Fig|2ji). In this case, we 
assume that /i jumps abruptly at y = between m for 
y < 0, and [i r at y > 0. The condition for the Majorana 
state to form is: 



/i|<B 2 -A 2 , fi 2 r >B 2 



(7) 



We match the wave function at y = 0, using the ansatz 
ip r oc e~ krV for y > and ipi oc e k ' y for y < 0. The 
Hamiltonian becomes: 

H = -(@{y)k r -Q(-y)ki)iuT z a z -pLT z + B(j x + AT x = 

(8) 



where ttfc Pj (i)± = A ± y B 2 — fJ- 2 ^ r an d the eigenvectors 
= e -(+)fe rfI) ±!/(i )e ±i9r(o i i,-ie ±ie -w) T /2 (9) 



with e % °rm = ^,, {l) /B + i^Jl - ^ 2 r{l) /B 2 . It is straightfor- 
ward to verify that ipr(i) ' ^ = (VV(i) ' are Majorana 
operators, with ip a simple c-number. Thus, we find that 
the wave function t^iv) of the Majorana state is 



2i sin 8 r ■ ip 



(OK 



)^l° )r + (e* 



y<0 

y>0 



(10) 

which exhausts all possibilities for isolated majorana 
states. 

Indeed we must note that when Eq is a spin gap, the 
gap Ei is due to pairing between spin-up electrons for 
positive p and spin-down electrons for negative p, remi- 
niscent of a one-dimensional p-wave superconductor [Io| . 
Recalling that vortices of a p-wave superconductor sup- 
port a zero-energy bound state 0, H, EH , we expect the 
formation of Majorana states when A changes sign (Fig. 
Hk). Due to the broken azimuthal symmetry, however, 
two inseparable Majorana states form where A vanishes. 

Next we discuss experimental realizations. The main 
requirement for our proposal to be feasible is a suffi- 
ciently strong spin-orbit interaction. Spin-orbit coupling 
in wires adiabatically connected to reservoirs was con- 
sidered long ago, both without electron-electron interac- 
tions [l2l ] and with interactions [l3T | in the framework of 
Luttinger-liquid theorv. Recently, this problem attracted 
renewed theoretical [3| and experimental interest, 
both with and without external magnetic field. 

Several candidate systems for quantum wires with 
spin-orbit interaction exist. In carbon nanotubes, spin- 
orbit coupling arises due to curvature effects [l6| . Here 
it is preferable to have a strong spin-orbit coupling along 
the direction of propagation, requiring that the tube is 
bent along its axis. Alternatively, one could introduce 
a strong electric field perpendicular to the axis. Per- 
haps a more promising candidate is a wire of InAs in the 
wurtzite structure which is known to have strong spin- 
orbit coupling \vi\ . The velocity u in the Hamiltonian Eq. 
(HJ) is related to the experimentally measured length scale 
Aso = lOOnm = mu and Ago = 250fiV = mu 2 /2 via 
u ~ /i2AsoAso ~ 7.6 x 10 6 cm/sec and m = h 2 /Xg Q 2A = 
0.015m e , with m e the free electron mass. Similar num- 
bers (with A = 280^ V) describe newly fabricated InSb 
wires, except with a large g-factor of ~ 50, compared to 
g ~ 8 in InAs, requiring only a small, relatively innocu- 
ous to the SC, magnetic field [l8| . 

The wire-Majorana states we envision, can be formed 
by spatial variations of the Zeeman field, the proximity- 
induced superconductivity, or, most importantly, the 
chemical potential, and will form near points where 
B 2 — (fi 2 + A 2 ) = 0. A varying chemical potential, as 
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in Fig [2Ji, for instance, can be achieved by gate elec- 
trodes capacitatively coupled to the wire. Tunneling ex- 
periments should provide the most direct signatures of 
the Majorana states [19j . 

Additional experimental signatures can be probed by 
controlling the phase of the pairing A in addition to the 
chemical potential. In particular, the configuration of 
Fig. [3] allows controlling the pairing phase on the left, 
center, and right sections independently; we denote these 
phases by <j)g, 4> c , <p r . The total Joephson current flowing 
between the three superconducting segments is rather in- 
tricate, and will be discussed in a separate publication. 
Since the Majoranas are localized when the distance be- 
tween them, L, is infinite the Jospehson current due to 
the Majoranas is zero. A straightforward first-order per- 
turbation analysis for finite L yields the energy splittings 
between the two Majorana states on the domain walls 
(c.f. Ref. d^). We find the Josephson energy associated 
with the Majorana fermions to be: 




— ) + E c cos(- 



0c). (11) 



= in the center region, and 
sides. Also, Ei r ~ E c « 

e -(B-A)L/u_ j n the s i m il ar 



E = Eg r cos( — 

here we assume that fi c 
fit = [x r = (i on the 

/ 2AB(A 2 -B 2 )(A 2 +[i 2 -B 2 ) 
V A 2 (A 2 + t i 2 -B 2 )-B t j. 2 (B+A) 

setup of the edges of a topological insulator 2jJ, l22| the 
E c , which is a result of the tunneling from the left and 
the right sections to the middle section is absent since 
the center region between the Majoranas is not in prox- 
imity to a superconductor. We notice two prominent 
effects. First, by letting <\> c = <f> T) for instance, we see 
that the Josephson current from the left superconductor 
is Att periodic. More interestingly, when we try to draw 
current from the center region, the current, proportional 
to E c sm((<j>e + </> r — 2<j> c )/2) is drawn equally from the 
left and right superconductors. Therefore this geome- 
try can serve as a Josephson transistor, since a change 
of (f>R determines part of the current between the left 
and the middle section. Related effects were discussed in 
Ref. (23 1 . 

Next, let us estimate the Josephson tunneling 
strengths in the setup of Fig. 02 with three independent 
superconducting substrates. If the substrates are prop- 
erly insulated from each other, only Josephson currents 
will be carried through the quantum wire via proximity; 
these will be given by Eq. (|I1[) , in addition to a contribu- 
tion proportional to A which is 2n periodic in <p r — <f> c and 
<f>e — <f> c . For the energy scales associated with the In As 
wires, we expect the critical current for the 2ir periodic 
portion to be of order AOnA, consistent with B,A^ IK. 
The 47r-periodic critical Majorana- Josephson currents, 
2eE c /H and 2eEi r /h, are a significant fraction of this 
number. For instance, for InAs parameters with fi c = 
and u r = M = 0.95 and A = 0.85, with B ~ IK, 
we obtain E c m 0.22Ke~ L ( B ~ A ^ u , corresponding to a 
maximum current eE c /h w 4nA ■ e - L / 3 A"« with L the 



FIG. 3. fi, B or A are tuned so that in the middle supercon- 
ductor we have a spin gap while the other superconductors 
have a superconduting gap. In this configuration, it is pos- 
sible to manipulate the two junctions separately by changing 
the superconducting phase difference between the neighboring 
regions. 



separation between the Majoranas. The unique flux pe- 
riodicity of the Majorana- Josephson currents can also be 
probed with low-frequency shot noise would also reveal 
the anomalous Josephson periodicity. 

In this manuscript we have shown that wires with 
strong spin-orbit coupling in proximity to a supercon- 
ductor host an interesting effective helical state. By tun- 
ing the superconducting gap A, the spin gap B, or the 
chemical potential fx, Majorana states can be created and 
detected in various experimental ways. By fabricating a 
set of gates over a network of such wires, we can imag- 
ine adiabatically creating Majorana pairs, moving, and 
even interchanging them along the network using pulse 
sequences in the gates. The non-Abelian character of the 
system should become apparent in such networks. Meth- 
ods to manipulate the Majorana modes and experimental 
consequences for the conductance will be the subject of 
a future manuscript. 

While finishing this manuscript we became aware of 
the preprint [24j which has some overlap with our results. 
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